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The current-induced magnetization dynamics of a spin valve are studied using a macrospin (single 
domain) approximation and numerical solutions of a generalized Landau-Lifshitz-Gilbert equation. 
For the purpose of quantitative comparison with experiment [Kiselev et al. Nature 425, 380 (2003)], 
we calculate the resistance and microwave power as a function of current and external field including 
the effects of anisotropies, damping, spin-transfer torque, thermal fluctuations, spin-pumping, and 
incomplete absorption of transverse spin current. While many features of experiment appear in 
the simulations, there are two significant discrepancies: the current dependence of the precession 
frequency and the presence/absence of a microwave quiet magnetic phase with a distinct magne- 
toresistance signature. Comparison is made with micromagnetic simulations designed to model the 
same experiment. 



I. INTRODUCTION 

In the ten years since Slonczevirski^ and Berger- dis- 
covered the quantum mechanical phenomenon of spin- 
transfer torque,'^ considerable evidence has accumulated 
that a spin-polarized current that passes through a thin 
ferromagnetic film can induce switching and/or preces- 
sion of the film's magnetization. Early experiments 
used multilayers,'* nanowires,^ small particle junctions,^ 
and point contacts? to infer the presence of this effect. 
The most convincing data came later from pillar-type 
"spin valves" with nanometer-scale transverse dimen- 
sions (Fig^ where a thin film non-magnet is sandwiched 
between two thin film ferromagnets^iS In the range of 
film thicknesses most commonly used, the magnetization 
M of the thick "fixed" layer and the magnetization m of 
the thin "free" layer lie in the plane of the film. Non- 
magnetic leads connect the spin valve to electron reser- 
voirs. 

Due to the phenomenon of giant magnetoresistance,^^ 
voltage measurements are sufficient to reveal that hys- 
teretic switching of m occurs as a function of the ap- 
plied current density J when a magnetic field H smaller 
than the coercive field is applied along the easy axis of 
the free layer. For larger values of H, it is believed 
that m exhibits one or more types of stable preces- 
sion as a function of J until the current density is large 
enough to induce switching. This conclusiorkiitiSiiiiiii^ 
is based on the experimental observation of narrow band 
microwave emission combined with calculations using a 
generalized Landau-Lifshitz-Gilbert (LLG) equation that 
predict precession of the free layer. Other observed dy- 
namical behavior includes telegraph noise that is inter- 
preted as rapid switching between two distinct states of 
magnetization>iSiiLi2iiS 

Several experimental groups have used the macrospin 
(single domain) approximation to propose "phase dia- 
grams" that identify the dynamical state of their spin 
valves as a function of J and i/ ,8^^17^^,20,21,22^ rj.^^^^ 



FIG. 1: Side view of a spin- valve (schematic). Anon-magnetic 
spacer layer is sandwiched between a "fixed" ferromagnetic 
film with uniform magnetization M and a "free" ferromag- 
netic film with uniform magnetization m. The leads on ei- 
ther side of the sandwich are non-magnetic. The free layer 
has thickness d and the direction of positive electric (negative 
electron) current J > is indicated. 



have also been purely theoretical studies of the LLG 
equation (generalized to include spin-transfer torque) 
using both macrospin niodels^iSSiSSil'^'^s. 29.30 g^j^j 
cromagnetics simulation a?^i?^i??i?'*i?^i?^i?Ti?8 that do not 
make the single-domain approximation. Unfortunately, 
it is difficult to extract a coherent picture from all this 
work because different authors make different choices for 
the physical effects they believe most affect the dynam- 
ics. There is not even unanimity amongst authors for the 
form of the spin-transfer torque itself. 

This state of affairs motivated us to perform a thor- 
ough study of the LLG dynamics of a model spin- valve for 
the purpose of a quantitative comparison with the data 
reported by Kiselev et alm^ for a Co/Cu/Co nanopillar. 
We make the macrospin approximation, but otherwise 
systematically examine the effects of different forms of 
spin-transfer torque, thermal fluctuations, spin-pumping, 
incomplete absorption of transverse spin current, and 
angle-dependent damping. We find that a "minimal" 
macrospin model can reproduce many (but not all) fea- 
tures of the experiment. The most important points of 
disagreement are the current dependence of the preces- 
sion frequency and the existence of a microwave quiet 
magnetic phase with a distinct magnetoresistance signa- 
ture. In light of these results, we comment on micromag- 




FIG. 2: The ellipsoidal cross section of the free layer (shaded) 
lies in y-z plane. We represent its magnetization by a 
macrospin m that can point in any direction. The fixed layer 
(not shown) is represented by a fixed macrospin M || z. 



netic simulationgS^i^ designed to model the identical set 
of experimental data. 

The plan of this paper is as follows. Section II de- 
scribes the macrospin models of interest and the gen- 
eralized Landau-Lifshitz-Gilbert equation we solve nu- 
merically. Section III presents results for a "minimal" 
model and compares them to the measurements reported 
in Ref. ilL. Section IV examines several variations of the 
minimal model within the context of the macrospin ap- 
proximation. Section V compares our results with micro- 
magnetic simulations. Section VI summarizes our results 
vis a vis experiment. An Appendix provides some details 
omitted from the main body of the paper. 



II. THE MACROSPIN MODEL 

Our macrospin model of the spin valve shown in Fig.Q] 
assumes that the magnetization is spatially uniform in 
both ferromagnetic layers with saturation value Mg. The 
fixed layer magnetization is M = M^z, but we allow the 
unit vector in the direction of the free layer magnetization 
m = m/Ms to point in any direction. In the coordinate 
system used here (Fig. |21), 
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(1) 



The experiments of interes^ii use a "free" ferromagnetic 
layer with a thickness cZ w 3 nm and an elliptical shape 
of dimensions about 130 nm x 70 nm. Under these con- 
ditions, magnetostatic shape anisotropy makes the y-z 
plane an easy plane for m. The z-axis is an easy axis 
in that plane. The control parameters are an external 
magnetic field H directed along ~\-z and an electric cur- 
rent J that is reckoned positive when negatively charged 
electrons flow from -I- a; to —x. 

We describe the dynamics of rh using a generalized 
Landau-Lifshitz-Gilbert (LLG) equatioufi*^ 



dm 
~dt 



-7m X [Heft + Ht] + am x 
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dm 



dt fJ,oM> 



■N. (2) 



It will be convenient to discuss each term in Eq. (0) in 
turn. 



A. Energy 

The first term on the right side of Eq. is a conven- 
tional magnetic torque with gyromagnetic ratio 7. This 
torque is driven by an effective field derived from the 
total energy E of the free layer with volume V: 



1 dE 



(3) 



Taking account of magnetostatics, the external field, and 
a uniaxial surface anisotropy, we show in Appendix I that 
E can be written in the form 



2E 



= hz cos" 6 + hy sin 9 sin (p + hx sin 9 cos 
— 2ft, cos 6*. 



(4) 

Here, h = H/Mg and the constants /ix, fty, and hz are 
computed in Appendix I using the free layer data given 
just below Eq. ^ and the material constants listed in 
Table I. 



B. Damping 

The "Gilbert damping" term am x rh in Eq. Q takes 
account of energy dissipation mechanisms such as cou- 
pling to lattice vibrations^ and spin-fiip scattering!^ 
The prefactor a is usually treated as a phenomenological 
constant (Table I) although it is not known whether this 
is a good approximation for situations where the ampli- 
tude of precessional motion is large. The Landau-Lifshitz 
approach to damping replaces the Gilbert term in Eq. Q 
by 



Am X (rh X Hcff). 



(5) 



The constant A can be calculated in some microscopic 
models?^ but a phenomenological treatment is almost 
universal. When N = in Eq. the Gilbert and 
Landau-Lifshitz expressions for the damping torque are 
known to be equivalent, at least formally.*'^ Section Hi PI 
gives a reason why we prefer the Gilbert form, but we 
performed calculations using both forms for purposes of 
comparison. No significant differences were found. 



C. Thermal Fluctuations 

The stochastic vector Ht in Eq. (0) is used to simulate 
the effect of finite temperature. Each Cartesian compo- 
nent is chosen at random from a normal distribution with 
a variance chosen so the system relaxes to a Boltzmaim 
distribution at equilibrium^^ Specifically, 

2kBTa 
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(6) 
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Quantity 


Values 


Ms 


0.127 X 10'^ A/m [llj 




1.6 T 


7(Co) 


2.4 X 10^ m/(A-s) [45] 


a{Co) 


0.01 [4^ 


Ku 


0.5 X 10"^ J/m^ 


s"/S(Cu) 


2.94 X 10^^ m"^ [48] 


V (Co/Cu) 


0.98 [48] 



TABLE L Quantity Values 



where «, j = x,y, z. We have confirmed numerically that 
this procedure does indeed produce a Boltzmann dis- 
tribution of energies at temperature T when N = in 
Eq. ^ . Details of our implementation of the stochastic 
contributions are indicated in Sec. IIII Al 



D. Spin- Transfer 

The quantity N in Eq. |5J) stands for one of several 
torque densities that arise from microscopic considera- 
tions of the transport of electrons through a spin valve. 
The most important of these is the spin-transfer torque 
density Ngti^*^ A variety of theoretical methods confirm 
the following picturei^ The current that flows through 
the spin valve shown in Figgis spin-polarized. Because 
M and m arc not collincar, the conduction electron spins 
that encounter the free layer generally possesses a com- 
ponent of angular momentum that is transverse to the 
magnetization of free layer itself. Realistic calculations 
show that this transverse component of angular momen- 
tum is largely absorbed by the ferromagnet , Since we 
describe the free layer as a uniformly magnetized parti- 
cle, the absorbed angular momentum generates a torque 
that appears on the right hand side of Eq. |j2J) . According 
to current theory^SSiSi^iiSiS 



2e a 



(7) 



where M = M/M^ and cos^ = m ■ M. 

The different forms of spin-transfer torque one finds 
in the literature correspond to different choices for r]{9). 
If one simply puts r]{9) ~ rjQ, the result is a "sine" ap- 
proximation to the torque because the remaining angular 
factors in Eq. (TJ give Ngt oc sin 6 (Fig. O . This form 
of the torque arises when there is spin-dependent scat- 
tering at the free layer interface and the polarization of 
the electron current that flows from the fixed layer to 
the free layer is independent of the orientation of the free 
layer. The prefactor 77(6*) is not constant if there is a dif- 
fusive component to the current anywhere and/or spin- 
dependent reflection occurs at the fixed-layer interface. 
To our knowledge, one or both of these effects is present 
in all transport theory calculations of Ngt- On the other 
hand, the corresponding sin^(0/2) approximation for the 
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FIG. 3: Various forms of (dimensionless) spin-transfer torque 
as a function of the angle 6 between the fixed layer and the 
free layer. The sine torque does not depend on the spin valve 
geometry. The symmetric Slonczewski (SS) and asymmet- 
ric Slonczewski (AS) torques are essentially identical for the 
standard spin valve geometry studied in this paper (AS/SS 
solid curve). The dashed curves show the difference between 
the symmetric and asymmetric Slonczewski torques for a ge- 
ometry discussed in Sec. llV Bl 



angular dependence of the magnetoresistance describes 
real spin valve data&S better than one would expect 
based on the transport theory predictions, to which we 
turn next. 

Building on his original workfi SlonczewskiS ap- 
plied an approximate form of magnetoelectronic circuit 
theory^^ to a spin valve with equal lead lengths and equal 
ferromagnetic layer thicknesses. He found 



A + Bcos9' 



(8) 



where q, A, and B are material and geometric parame- 
ters. We will call this the symmetric Slonczewski (SS) 
approximation for the torque. For a general spin valve 
geometry, it turns out thali^^i^ 



q- 



A + Bcos9 A - B cos 9 



(9) 



The present authors have shown that Eq. © gives 
quantitative agreement with calculations of the spin- 
transfer torque based on the Boltzmann transport equa- 
tion for a wide variety of spin valve geometries. "'^^ We 
will call this the asymmetric Slonczewski (AS) approx- 
imation. One of the solid curves in Fig. O shows that 
the symmetric and asymmetric Slonczewski torques are 
essentially identical for the particular spin valve geome- 
try we use to model the experimental sample of Ref. 
(see Sec. IIII|l . The two dashed curves show the differ- 
ence between the the symmetric and asymmetric Slon- 
czewski torques for a spin valve geometry we will discuss 
in SecHvB 

Spin-transfer torque accounts for non-equilibrium pro- 
cesses that cannot be described by an energy functional. 
This means that Ngt does not produce an effective field 



3 



like Eq. Q and no damping of spin-transfer dynamics oc- 
curs if Heff = and the Landau-Lifshitz form Eq. ^ is 
used for damping. On the other hand, if one beheves that 
it must be possible to influence spin-transfer driven mo- 
tion by transferring energy to other degrees of freedom, 
it is necessary to the use the Gilbert form of damping in 
the magnetization equation of motion. This is what we 
do in Eq. ||2J). 



with 



hq^^ dm 
AnS dt 



(13) 



S is the cross-sectional area of the free layer. Table I gives 
numerical values for the parameters v and g^^ (defined 
in Ref. l68j) for the Co/Cu/Co spin valve of interest to us 
here. 



E. Current-Induced Effective Field 

First principles calculations^Si^ show that the absorp- 
tion of a transverse spin current at a ferromagnetic in- 
terface is not 100% efficient. Part of the fraction that 
survives gives a small correction to ri{d) in Eq. Q. The 
remainder is polarized perpendicular to both m and M 
and contributes a torque density on the free layer of the 
form 

Ncff-?7(0)/3:^^mxM. (10) 
2e d 

In the circuit theory language of Brataas et al.^^ this term 
is described by the imaginary part of the mixing conduc- 
tance. Evidently, Nes produces motion of m identical 
to that produced by an effective external field oriented 
along the magnetization direction M of the fixed layer. 
This contribution is usually neglected because the cited 
calculations find (3 « 0.05. We include it here because at 
least one experimenliS^ has been interpreted as demon- 
strating that [3 « 0.20. 



III. MINIMAL MODEL 

This section compares LLG simulation results with the 
experimental results reported in Ref. Illl Our "minimal" 
model is Eq. Q with N = Ngt from Eq. (TJl and the 
asymmetric Slonczewski (AS) choice in Eq. lO for rj{6). 
This model takes account of magnetostatic and surface 
anisotropy, an external magnetic field, current-induced 
spin-transfer torque, Gilbert damping, and thermal fluc- 
tuations. Most of our calculations use a spin valve geom- 
etry (see Fig.^ designed to mimic the nanopillar samples 
studied by Kiselev et ali^: 

Cu(80 nm)/Co(40 nm)/Cu(thin)/Co(3 nm)/Cu(10 nm). 

The notation Cu(thin) indicates that the thickness of the 
spacer layer is immaterial as long as it is smaller than 
the mean free path in copper. The precise choice of lead 
lengths is subject to uncertainty due to the approxima- 
tions needed to model finite width and reservoir effects 
in a one-dimensional Boltzmann equation calculation of 
spin valve transport^ 



F. Spin Pumping 



A. Computational Details 



A final contribution to the torque on the free layer 
comes from a phenomenon called "spin pumping" . Since 
a spin polarized current incident from a non-magnet can 
produce magnetization dynamics in an adjacent ferro- 
magnet, it is not unreasonable that motion of the mag- 
netization of a ferromagnet can influence the spin current 
in an adjacent non-magnet. The most prominent effect 
is the injection of a spin current into the non-magnet 
whenever the magnetization moves. One consequence of 
the injected spin current is a back-reaction torque that in- 
creases the damping of the spin motioniSi*SiS This effect 
has been confirmed by experimentsi^^i2S*2SiSi The torque 
density due to spin-pumping is given by Tserkovnyak et 
alM as 



where 



TCXCh 



M — j/m cos 6 
1 — v"^ cos^ 



(11) 



(12) 



We solve the stochastic LLG equation using the 
Ito calculus^^ and a numerical method described by 
Milshtein.^". The simulations proceed by fixing the exter- 
nal field H and sweeping the current density J in steps 
of size 5J . Before changing to the next value of J, we 
integrate the LLG equation for a "waiting time" t* using 
N time steps of length 5t. After each time step, we use 
the instantaneous value of the angle 6 between M and m 
and the results of Ref. to evaluate the instantaneous 
magnetoresistance R{9). A time-average over these N 
values gives the resistance we report for each J . 

Fig. 0] shows the calculated high-field and low-field 
magnetoresistance as a function of J for three values of 
the simulated sweep rate SR = 6J/t*. The curves in this 
figure are averages over 20 realizations of the stochas- 
tic simulation. In each realization, the system switches 
abruptly at a particular value of current between states 
with distinctly different magnetoresistance (to be dis- 
cussed below). Since the switching current depends on 
the realization, an average over essentially vertical tran- 
sitions at slightly different switching currents gives the 
not-quite-vertical lines seen in the figure. As expected. 
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7(10^ A/cm2) 



FIG. 4: High field (upper panel) and low field (lower panel) 
magnetoresistance as a function of current density sweep rate 
(SR) in units of A/(cm^-s). The up arrows identity the parts 
of the hysteresis loops traced out when J is scanned from neg- 
ative values to positive values. The down arrows correspond 
to scanning from positive values to negative values of J. 



the hysteresis loops close as the sweep rate decreases. 
Less obviously, the rate of closing is much greater at 
high field than at low field. It is important to appre- 
ciate that the slowest sweep rate we can practically use 
in our simulations [10^^A/(cm^-s)] is still five orders of 
magnitude faster than the sweep rate used in the Cornell 
experimentsiS 

For fixed values of H and J, the N values of resistance 
collected between t — and t = t* constitute a time 
series for the resistance. Spin valves are Ohmic devices, 
so the Fourier transform of this series is proportional to 
the associated power spectrum. We use this numerical 
data below to compare with the microwave noise data 
reported in Ref. ITlL 



B. J-H Phase Diagrams 

Fig. compares spin valve "phase diagrams" at T = 
3 K, T = 300 K, and T = 3000 K for our minimal model. 
The diagrams were constructed by sweeping the current 
twice (once increasing the current and once decreasing 
the current) for each value of H. There is some noise 
at higher temperature because we did not average over 
multiple realizations of the simulation. Solid lines divide 
each diagram into phase fields with labels like A, B and 
A/B. The latter means that the field is occupied by phase 
A when the current is scanned from left-to-right in the 
diagram and by phase B when the current is scanned 
from right-to-left. Thus, a label like A/B is a signal that 
hysteresis is present. 

The phase fields in Fig. [51 are labelled P (parallel), 
AP (anti-parallel), IPP (in- plane precession) and OPP 
(out-of-plane precession). The static P and AP states 
are labelled by the relative orientation of m and M. 
The precessing states are identified from the microwave 



power (not shown) as described above. IPP denotes a dy- 
namic state where m precesses symmetrically (or nearly 
so) around an axis that lies in the y-z easy plane. OPP 
denotes a dynamic state where m precesses symmetri- 
cally (or nearly so) around an axis that does not lie in 
the easy plane. Sec. IIII Dl describes these states in more 
detail. 

We focus first on the 3 K diagram. This is similar 
(but not identical) to T = K diagrams published by 
othersiiSi using the symmetric Slonczewski spin-transfer 
torque. Using sharp peaks in the measured noise power 
spectrum to identify states of stable precession, Kiselev 
et al. pointed out the topological similarity between their 
computed T — K phase diagram and their measured 
T — 300 phase diagramiii 

When H exceeds the coercive field, our 3 K phase di- 
agram shows hysteresis for the P^IPP and OPP<-^AP 
phase transitions. The experiment shows no hystere- 
sis in this regime (see below). At low field, the P — > 
AP transition occurs abruptly while the reverse-current 
AP — > P transition does not. Instead, there is a long, 
skinny, triangular-shaped P /IPP phase field within which 
the magnetization m exhibits stable, elliptical preces- 
sion around the — z axis. The precession amplitude in- 
creases as the current becomes more negative. The sys- 
tem crosses the phase boundary into the P phase when 
the precession angle between m and — z exceeds 90° and 
the vector m spirals irreversibly toward z. We will see 
below that this asymmetry has its origin in the details of 
the dependence of the spin-transfer torque on the angle 
between the free layer and the fixed layer. 

We draw special attention to the lower limit of the 
OPP/AP phase field in the 3 K phase diagram. The 
perfectly horizontal portion of this phase boundary is an 
artifact of the current scanning mode used to generate the 
diagram. If we fix J and scan the external field H from 
large values to small values, the OPP phase does not give 
way to the AP phase until the dashed line in the diagram 
is crossed. The exact shape of this boundary depends on 
the i/-scan rate, but it is reasonable to suppose that 
the corresponding horizontal phase boundary found in 
Ref. may also be an artifact of the method of taking 
data. 

Reading Fig. [S] from left to right shows that the re- 
gions of hysteresis shrink as the temperature increases. 
More details can be seen in the line scans of Fig. 
The effect of increasing temperature is very similar to 
the effect of decreasing the current density sweep rate in 
Fig. ^ Indeed, since our simulated current sweep rate 
is always much faster than experiment, the temperatures 
indicated on the phase diagrams in Figure [S] must be re- 
garded as nominal. The true phase diagram at each tem- 
perature we show would exhibit less hysteresis. Equiva- 
lently, each panel actually corresponds to a lower physical 
temperature than the temperature we quote. Thus, our 
T = 3000 K diagram indicates (qualitatively) how the 
300 K phase diagram might look if we could use cur- 
rent sweep rates comparable to those used experimen- 
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FIG. 5: Minimal model dynamic phase diagrams for a Cu(80 nm)/Co(40 nm)/Cu(thin)/Co(3 nm)/Cu(10 nm) spin valve. Left 
panel: T = 3 K; middle panel: T = 300 K; Right panel: T = 3000 K. For fixed H, a bistable region labelled A/B exhibits the 
A state when J is scanned from left to right and the B state when J is scanned from right to left. The correspondence needed 
to compare with Ref . is 10* A/cm^ « 10 mA. The dashed curve is the OPP-^AP phase boundary for a field scan from 
large H to small H at fixed J. 
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FIG. 7: Telegraph noise in a time-series of the magnetore- 
sistance. The data was collected on the high-temperature 
boundary between the AP phase and the OPP phase {fJ-oH = 
0.04T, J = 0.7 X 10*A/cm^). 



FIG. 6: Temperature dependence of the magnetoresistance 
obtained by averaging the results of 20 realizations of the 
stochastic simulation. Upper panel: fioH = 0.05 T; Lower 
panel: fioH = 0.002 T. 



tally. We note also that substantial Joule heating occurs 
in real spin valve samples, perhaps 15 K to 20 K per 10^ 
A/cm2,ia 

Our highest temperature simulation shows P«-^AP hys- 
teresis when H is small and complete reversibility when 
H is large. This resolves the disagreement between 
theory and experiment noted above. Moreover, state- 
to-state switching characterizes every reversible phase 
boundary. Fig. [7| illustrates this for the dynamics of 
switching between the anti-parallel AP state and out- 
of-plane precession (OPP). The time-series data for the 
magnetoresistance was collected on the OPP/AP phase 
boundary at 3000 K. Clearly, the system switches back 



and forth between the AP state (small fluctuations 
around unit normalized resistance) and the OPP state 
(periodic oscillations of the normalized resistance be- 
tween zero and one). Experiments show precisely this 
sort of telegraph noiseiSii^iiSiiS if we replace the full scale 
excursions of the OPP resistance with small scale fluctu- 
ations around the average resistance of the OPP state. 
We are not aware of experiments that study the tele- 
graph noise at our 300 K IPP/OPP boundary or our 
3000 K P/OPP boundary. Indeed, at the latter, our sim- 
ulations actually show random switching between three 
states: AP,IPP, and OPP. 

The variations of the computed resistance near the 
3000 K P/OPP and OPP/AP phase boundaries lead 
to two peaks in the differential resistance, dV/dl — 
R + I dR/dl . These agree well with the peaks in dV/dl 
observed experimentally. On the other hand, Kiselev et 
al.^^' identify a "W"-phase that is completely absent from 
our 3000 K phase diagram. However, the experimen- 
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tal W-phase field appears exactly where our model pre- 
dicts OPP/AP phase bistability at 3 K and 300 K (two 
left panels of Fig. . The experimental W-phase is mi- 
crowave quiet above the experimental low frequency cut- 
off (0.1 gHz) and it exhibits a magnetoresistance that 
is slightly, but distinctly, smaller than that of the AP 
configuration. This would occur in our macrospin model 
if the free layer were frozen into a static configuration 
with m neither parallel nor anti-parallel to M. We will 
return to the W-phase when we discuss micromagnetic 
simulations in Section Ivl 

Quantitatively, our calculated coercive field is about 
half the experimental value. This discrepancy may re- 
flect an inaccurate description of the shape (and therefore 
the magnetostatic anisotropics) of the free layer. Another 
contributory effect is our complete neglect of dipolar cou- 
pling to the fixed layer. At low T, we also find that the 
magnitude of the critical current J+ for the P^AP tran- 
sition is much greater than the magnitude of the critical 
current J~ for the AP— >P transition. Experiments show 
that and J~ are more symmetric around zero current. 
The reason for our calculated behavior is"'^ the difference 
between 77(0) and ?7(180°) for the Slonczewski torques in 
Fig. O Our simulations at 3000 K more nearly resemble 
the experiments at 300 K because J+ decreases strongly 
with temperature while J~ is nearly temperature inde- 
pendent. 

C. Precession Frequency 

The gray/color scale in Fig. |S1 quantifies the relative 
microwave power (on a logarithmic scale) at frequency / 
as a function of magnetic field for two values of current 
density J. The numerical data was obtained by Fourier 
transforming our simulated time series data for the mag- 
netoresistance. The narrow bands of peak microwave 
power trace out the frequency uj{H) of stable precession 
(and its harmonics). The left panel corresponds to small- 
amplitude, noise-driven, in-plane precession at a value of 
current just before the parallel phase becomes unstable 
to steady in-plane-precession. As with the experimen- 
tal datafiiii^ ^{H) in this regime can be described by 
the Kittel equation (black curve) for thin film magnetic 
resonance!^ In the notation of Appendix I, the resonance 
frequency is 

iOK = 1^J[H + {hy - hz)Ms][H + (hx - hz)M,]. (14) 

Our simulation data agree with 2ll!k because the peri- 
odicity of the resistance is twice the periodicity of the 
magnetization oscillation frequency. No analytic theory 
is available for comparison with our results at a higher 
value of J (right panel) where large amplitude out-of- 
plane precession occurs. But our results do show the 
same relative magnitude and if-dependence as seen in 
the experiments. 

Fig.|51shows the relative microwave power at frequency 
/ as a function of increasing J (top panel) and decreas- 




FIG. 8: Relative microwave power at different frequencies as 
a function of field H at T — 3 K. The gray/color scale is 
logarithmic. Left panel: J — 0.3 x 10*A/cm^; the black curve 
is twice the resonance frequency given by the Kittel equation 
Eq. lO. Right panel: J = 1.0 x 10*A/cm^ The muhiple 
curves are the fundamental and its harmonics. 



ing J (bottom panel). Similar plots for comparison with 
experiment have been presented by others using the sym- 
metric Slonczewski torque^^ and the sine torque.'^" The 
zero power regions at low and high J correspond to the 
static P and AP magnetization states. In between, the 
narrow bands of peak microwave power trace out a;(J) 
(and its harmonics) for stable precession. Just above the 
limit of the parallel state, there is a very narrow range of 
in-plane precession where tt'(J) decreases monotonically. 
At slightly higher J, the system evolves to a state of out- 
of-plane precession (OPP) where uj{J) first increases and 
then decreases. Comparison of the two panels in FigureEl 
illustrates the hysteresis present at this low temperature. 

Our results for w( J) do not agree with observations for 
real spin valves. Putting aside the fact that no hys- 
teresis is seen in the experiments (which we attribute to 
the current sweep rate as discussed above), the experi- 
mental data always show that oj decreases as J increases. 
Naively, it is as if in-plane precession persisted all the way 
to the anti-parallel state with no intervening state of out- 
of-plane precession. This is a serious issue because, in our 
model, in-plane precession occupies an extremely small 
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FIG. 9: Microwave power at different frequencies as a function 
of J at T = 3 K and fioH = 0.05 T. The gray/color scale is 
logarithmic. Left panel: current scan from left to right. Right 
panel: current scan from right to left. The narrow bands of 
peak power represent the frequency i^{J) of stable precession. 
The black traces are the magnetoresistance in arbitrary units. 



portion of the J-H phase diagram. 



D. Precession Trajectories 

To help shed hght on our simulation results for oj{J), 
it is instructive to analyze the relationship between this 
quantity and the trajectory of the tip of m on the unit 
sphere. We will call this the orbit of the precessional 
motion. Without loss of generality, we set hy = hz — 
and retain only the external field H ~ Hz and the 
local, out-of-plane demagnetization field = —hxfhxX. 
It is crucial that the magnitude Hd changes along the 
trajectory because the component rhx of m changes along 
the trajectory. 

It is common to think of precession as the steady mo- 
tion of a vector on a cone that makes a small angle with 
respect to its symmetry axis. The orbit in this case is 
a circle. The precessional states in the present problem 
are more complicated. Fig. IIUI shows two large ampli- 
tude, saddle-shaped, in-plane precession (IPP) orbits for 




FIG. 10: In-plane precession (IPP) orbits for two nearby val- 
ues of J. The thick segments are points on each orbit where 
the magnitude of the demagnetization field — —hxthxit 
is smaller than the magnitude of the external field H — Hz. 
The y-z easy plane and the equatorial circle of the unit sphere 
in the x-y plane are indicated as guides to the eye. 



AS/SS 




S Sine 

FIG. 11: Out-of-plane precession (OPP) orbits. The geom- 
etry is the same as Fig. 1101 except that the north (N) and 
south (S) poles of the unit sphere are indicated. The orbit 
labelled "AS/SS" is produced by Slonczewski's spin-transfer 
torque. Along the thick segment, the spin-transfer torque and 
demagnetization field torque point in opposite directions. The 
orbit labelled "sine" is produced by a sine-type spin-transfer 
torque for the same value of J. 



two nearby values of J. We call these "in-plane" preces- 
sion modes because each orbit moves symmetrically (or 
nearly so) around an axis (the 2:-axis) that lies in the easy 
y-z plane. 

Let us partition each orbit into two segments. The 
short thick segments lie near the easy plane where the de- 
magnetization field Hd is smaller than the external field 
H. Along these segments, the orbital azimuthal angle 
(j) precesses mainly around H with angular speed jH. 
Along the remaining segment of each orbit, Hd is larger 
than H and the orbital polar angle 9 precesses mainly 



8 



around with angular speed ^Hd- Fig. 1101 shows that 
the angular range swept out by both the thick and thin 
segments increases as the current density (and the spin- 
transfer torque) increases, i.e., the total arc length of the 
orbit increases. Since the orbital speeds change very lit- 
tle with J, we conclude that the orbital period increases 
as current density increases. This implies that lij{J) is a 
decreasing function for in-plane-precession orbits. 

As J continues to increase, the apices of the two thick 
segments of the in-plane precession saddle orbit approach 
and then touch one another near the negative z-axis. 
When this occurs the orbit bifurcates into two elliptical 
orbits, each centered on an out-of-plane axis not far from 
the cc-axiSfiiii^^iSi Precessional states at higher current 
density correspond to one or the other of these out-of- 
planc (OPP) trajectories, e.g., the AS/SS orbit in Fig.lTTl 
This orbit precesses mostly around H^. Spin-transfer 
torque tends to push the orbit away from the easy plane 
in the northern unit hemisphere. The effect on the orbit 
in the southern unit hemisphere is more complex. The 
net result is that the "center" of the orbit moves away 
from the easy plane. In other words, as the current den- 
sity increases, the component "nix of m increases, which 
increases Hd, and thus increases the frequency ^Hd of 
the orbit. 

Along the thick segment of the out-of-plane orbit, the 
spin-transfer torque and the torque from the demagne- 
tization field point in (nearly) opposite directions. This 
means that the net torque, and thus the orbital speed 
along that segment, decreases as J increases. Eventually, 
this slowing down overwhelms the speeding up described 
just above and the total orbital period begins to increase. 
This is why u}{J) decreases for the largest values of J 
where precession occurs in the top panel of Fig. 1^1 



IV. BEYOND THE MINIMAL MODEL 

There are two major discrepancies between the Cornell 
experiments^ and our minimal model results: the varia- 
tion of the precession frequency w with current density J 
and the absence of a microwave quiet "W-phase" . Within 
the context of the macrospin model, we examined several 
variations of our model, mostly with the hope they would 
improve the agreement between theory and experiment. 
We studied the influence of (A) a sine-type spin trans- 
fer torque for the standard geometry; (B) asymmetric vs. 
symmetric Slonczewski spin-transfer torque for a special 
asymmetric geometry; (C) the current-induced effective 
field that arises due to incomplete absorption of trans- 
verse spin currents; (D) spin-pumping; and (E) angle- 
dependent Gilbert damping. 

A. Sine Spin- Transfer Torque 

Fig. shows the geometry-independent "sine" torque 
that is widely used in the literature. The top left panel 



in Fig. ^1 shows the 3 K phase diagram when this sine 
torque replaces the AS/SS torque. Several differences 
with the corresponding Slonczewski torque phase dia- 
gram (left panel of Fig. [SJ should be noted. 

First, with a sine torque, the low-field P^AP transi- 
tion is mediated by in-plane precession in the same way 
that precession mediates the AP^P transition for both 
the sine and Slonczewski torques. This occurs because 
sin 9 is symmetric around 9 — i: /2 while the minimal 
model torque is not. Second, the sine torque generates 
no hysteresis in the high-field transitions P<-^IPP and 
IPP^OPP. Third, the lower limit of the OPP phase 
boundary determined by a field scan from large H to 
small H (dashed curve) greatly reduces the size of the 
AP phase field compared to the Slonczewski case. This 
feature does not appear to have been noticed in previous 
discussions of this phase diagrami ^^i^^ 

Unlike the Slonczewski torque, increasing current or 
field eventually drives the sine torque model to a tran- 
sition from out-of-plane precession (OPP) to a static 
phase where the macrospin m is "fixed" (F) at some an- 
gle between zero and tt. This is intriguing because the 
magnetoresistance and microwave power characteristics 
of this phase match exactly to those of the experimen- 
tally observed "W-phase". Unfortunately, the location 
of the F phase in the sine torque phase diagram does 
not agree with the location of the W-phase in the experi- 
mental phase diagram (see the penultimate paragraph of 

sec Hira 

The bottom left panel in Fig. [H shows the 3000 K 
phase diagram for the sine torque macrospin model. 
Compared to the corresponding diagram for the minimal 
model (right panel of Fig. [SJ, thermal effects eliminate 
OPP/AP bistability only when J is small. Hysteresis 
between these phases remains when J is large. In this 
sense, the sine torque model is more resistant to thermal 
fluctuations than the Slonczewski torque model. 

Finally, the out-of-plane precession frequency for the 
sine torque is a strictly increasing function of J. The 
argument is similar to the one in Sec. IIII Dl for the Slon- 
czewski torque. However, as the current increases, the 
different angular dependence of the sine torque causes 
the orbit of m to push steadily away from the easy plane 
everywhere and contract on the unit sphere (see Fig.lll|l. 
The frequency ^Hd increases monotonically because the 
demagnetization fleld Hd increases. The outward motion 
and areal contraction of the orbit continues as the current 
increases until the orbit area shrinks to a single point on 
the unit sphere. This is the signature of the flxed (F) 
phase. 



B. AS vs. SS Spin- Transfer Torque 

The AS/SS curve in Fig. |21 shows that the asymmetric 
Slonczewski torque used in our minimal model is essen- 
tially identical to the symmetric Slonczewski torque for 
the spin valve geometry of Ref. 11. This is not always the 
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Jq(10^ Aycm2) Jq{10^ A/cm2) ( 1 0^^ A/cm2) 

FIG. 12: Phase Diagrams at 3 K (top panels) and 3000 K (bottom panels). Left panel: sine torque; middle panel: sine torque 
plus current-induced effective field; right panel: Slonczewski torque plus current-induced effective field. For fixed H, a bistable 
region labelled A/B exhibits the A state when J is scanned from left to right and the B state when J is scanned from right to 
left. The dashed curves are the OPP— >AP phase boundaries for a field scan from large H to small H at fixed J. 



case. For example, compared to the geometry of Kiselev 
et al.,^^ a spin valve with film thicknesses, 

Cu(10 nm)/Co(40 nm)/Cu(thin)/Co(3 nm)/Cu(180 nm), 

is very asymmetric. The Cu/Co bilayers on opposite 
sides of the spacer layer are very different: the left bi- 
layer is mostly ferromagnet, the right bilayer is mostly 
non-magnet. The difference between the symmetric and 
asymmetric Slonczewski torques for this geometry is still 
small (compare the two dashed curves in Fig. (SJ. Nev- 
ertheless, it is large enough to produce small-angle, in- 
plane precession that "rounds" the low-to-high resistance 
jump in the AS hysteresis curve at the P— >AP transition 
in the lower panel of Fig. ^1 The same panel shows sim- 
ilar precessional rounding for the sine-type spin-transfer 
torque. However, the latter rounding disappears when 
shape anisotropy is turned off. The corresponding round- 
ing for the asymmetric Slonczewski torque does not dis- 
appear when shape anisotropy is turned off.^** 

Fig. ^] also shows that the critical current J+ for 
the P— *AP transition differs for all three spin-transfer 
torques while the critical current J^T for the reverse 
AP— >P transition distinguishes only the sine torque. 
This is a consequence of the fact that J+ ( J~) is inversely 
proportional the slope of the torque function r](9) sin 9 
(plotted in Fig.EJ at 6* = (6* = tt)} The fact that J+ is 
smaller for the AS torque than for the SS torque suggests 



\IqH = 0.05 T 




0.0 „ T 0.5 

7(10* A/cm2) 



FIG. 13: Spin valve hysteresis at 3 K for three choices of spin 
transfer torque using a geometry (defined in Sec. llV 1^ chosen 
to emphasize the difference between the AS and SS torques. 
Lower panel: H = 0.03T; Upper panel: H = 0.05T. 



that an asymmetric geometry like the one above may be 
desirable for some applications. 
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Current-Induced Effective Fields 



We mentioned in Sec. Ill El that a current-induced 
torque that acts hke an effective external magnetic field 
can arise due to incomplete absorption of a transversely 
polarized spin current at the interface between the spacer 
and the free layer. A recent experimental report^^ has 
been interpreted by its authors to mean that the size of 
this torque — displayed in Eq. ifTTijl — is much larger than 
theoretical estimates. Accordingly, Fig.ll2lshows how the 
spin valve 3 K phase diagram changes if we augment the 
sine torque (top middle panel) and the minimal model 
torque (top right panel) by an effective field torque that 
is 20% of the spin-transfer torque, as suggested by this 
experiment. For comparison, the bottom middle and bot- 
tom right panels respectively show the phase diagrams for 
these two situations at 3000 K. We find that the topology 
of the phase diagram does not change, although the pre- 
cise positions of the phase boundaries do. The qualitative 
behavior of the precession frequencies is not affected. 



D. Spin-Pumping 

The original discussion of spin-pumping^^ focused on 
the enhancement of Gilbert damping that occurs when 
a normal metal is in intimate contact with a precessing 
thin film ferromagnet. Subsequent workM has empha- 
sized that the torque due to spin-pumping is generally 
of the same order of magnitude as spin-transfer torques. 
Since the analytic form of the torque Ngp in Eq. Ullfl 
differs considerably from simple damping for, e.g., large- 
angle, out-of-plane precessional motion, this raises the 
possibility that spin-pumping alters the dynamical be- 
havior of a spin valve more profoundly than merely en- 
hancing the Gilbert damping. 

In the small angle limit, the parameters for intrinsic 
Gilbert damping and spin-pumping suggested in Ref. ^3 
for the Co/Cu/Co system (Table I) produce a total ef- 
fective Gilbert damping of = 0.148. Therefore, 
in Fig. El we compare the phase diagram at 3 K ob- 
tained with our minimal model (no spin pumping) us- 
ing a = Qfoff (solid lines) with the phase diagram ob- 
tained including both (reduced) Gilbert damping and 
the spin-pumping torque density Ngp (dashed lines) . The 
small differences we find between the two show that spin- 
pumping does not much affect the sort of precessional 
motion produced by our minimal model. There is also 
no qualitative change in the current dependence of the 
precession frequency. We conclude that neither effective 
fields nor spin pumping affects improves the agreement 
between experiment and the minimal model. 
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FIG. 14: Minimal model phase diagram at 3 K with Gilbert 
damping only (solid lines) and spin-pumping with reduced 
Gilbert damping (dashed lines). See text for discussion. 



nance or Brillouin light scattering experiments p£i Since 
the magnetization is never tilted far from equilibrium 
in these situations, it is relevant that Back et al."^^ re- 
ported an effective increase in the damping constant for 
large magnetization rotation angles in cobalt films under 
pulsed field conditions. More recent pulsed-field experi- 
ments have been successfully analyzed using conventional 
Gilbert damping^S*2i Since transient magnetization dy- 
namics is a common feature of all these experiments, it 
is not obvious that the results address the validity of 
the constant a approximation for large-angle, steady pre- 
cession of the sort discussed in this paper. With one 
exceptioufSi we are unaware of any models that allow the 
damping constant to vary during the course of precession. 

A glance back at Fig. |51 shows that the precession fre- 
quency u>{J) decreases monotonically (as seen in exper- 
iment) for in-plane-precession and also for out-of-plane 
precession near the boundary with the anti-parallel state. 
Using the damping "constant" ansatz 



(15) 



we have been able to produce dynamical phase diagrams 
with (i) only in-plane-precession or (ii) only out-of-plane 
precession with (in both cases) the "correct" behavior for 
ti'(J), or nearly so. Unfortunately, other features appear 
in the simulated phase diagram that do not appear in 
the data, e.g., an extended region of precession when the 
external field H is less than the coercive field. It may 
be possible to fine-tune the phase diagram to the desired 
form, but this seems unwarranted without some justifi- 
cation for Eq. (|15|l and the relative paucity of theoretical 
information about damping far from equilibrium^ 



E. Angle-Dependent Gilbert Damping 



MICROMAGNETICS 



The numerical value of the Gilbert damping constant 
a in Eq. [3 is usually extracted from ferromagnetic reso- 



Our macrospin approach to the Landau-Lifshitz- 
Gilbert (LLG) equation replaces the free layer magne- 
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tization m(r) by a constant vector m. The numerical 
method of micromagnetics^Si is a better approximation 
to reahty because it retains spatial gradients of m(r) 
down to a fixed minimum length scale. Treating the spa- 
tial variation of the magnetization allows the inclusion of 
two effects that we have neglected because they do not 
contribute when the magnetization is uniform. First, the 
experimental samples are polycrystalline so we ignore any 
intrinsic magnetocrystalline anisotropy because the non- 
uniform effective fields tend to average to zero over the 
whole sample. Second, the Oersted magnetic field pro- 
duced by the current itself largely averages to zero over 
the whole sample. While these two effects are not im- 
portant for the macrospin dynamics we have considered 
here, full simulations^"^ show that the inhomogeneities in 
the magnetization that result can be quite important. 

There is an emerging consensus amongst micromag- 
netic practitioner o?^i?^i??i?'^i?'^i??i?Ti?? that the inclusion 
of spin-transfer torque excites incoherent spin waves in 
the free layer (and thus inhomogeneous magnetization) if 
the current density is sufficiently great to induce switch- 
ing and microwave emission. However, since the method 
takes full account of local exchange and non-local mag- 
netostatics, systematic survey calculations of the sort we 
have presented in this paper are prohibitively expensive, 
even for systems as small as a spin valve free layer. Us- 
ing a torque density Ngt oc sm9 to model spin-transfer, 
Ref. 1^ reports finite-temperature micromagnetic simu- 
lations designed to mimic the experimental conditions 
reported by Kiselev et ali^ Like our Fig. [3 the cal- 
culated phase diagram agrees topologically and semi- 
quantitatively with the experimental diagram, with one 
important improvement. The micromagnetics simulation 
identifies the experimental "W-phase" with a dynamical 
phase field where vortices of magnetization continuously 
form and annihilate. The calculated noise power in this 
regime is concentrated at very low frequency and thus 
appears to be microwave quiet in the experiment. 

It seems plausible that the extra degrees of free- 
dom present in the micromagnetic approximation allow 
the bistable OPP/AP phase present in our macrospin 
model to break up into a spatially inhomogeneous state. 
The same logic suggests that the fixed F phase of the 
macrospin sine torque phase diagram breaks up similarly 
into a state of inhomogeneous magnetization. On the 
other hand, the micromagnetic simulation produces a 
nearly horizontal line for the AP phase boundary that 
we identify as an artifact of the current-scanning mode 
of data collection. 

It is worth noting that the microwave power output 
in this micromagnetic calculation is nearly independent 
of J in the high-field precession regime (as we find for a 
macrospin) while the power output observed experimen- 
tally varies considerably in this part of the phase dia- 
gram. The precession frequency w( J) in the same regime 
has been studied by Berkov and GornS, who find that 
they are able to qualitatively reproduce the experimen- 
tal frequency dependence with a highly inhomogeneous 



magnetic state. The inhomogeneous state arises from the 
non-uniform Oersted field and an (assumed) random dis- 
tribution of granular magnetocrystalline anisotropy. It 
would be interesting to discover if this type of micro- 
magnetic simulation can produce resonance linewidths 
with Q « 100 as observed in the most recent nanopillar 
experiments^ 



VI. SUMMARY & CONCLUSION 

We have studied the Landau-Lifshitz-Gilbert (LLG) 
dynamics of a single macrospin as a model for current- 
driven magnetization motion in the free layer of a spin 
valve. We parameterized our model specifically to com- 
pare our results with those reported by Kiselev et al}-^ 
for a Cu/Co/Cu nanopillar. Due to the simplicity of the 
model (compared to micromagnetics), we were able to 
explore systematically the effects of temperature, spin- 
pumping, current-induced effective fields, various forms 
of spin-transfer torque, and angle-dependent damping. 
We focused most of our attention on a "minimal model" 
where Slonczewski's spin-transfer torque supplements the 
terms usually found in the LLG equation. 

Low-Field Behavior: Our minimal macrospin model 
captures the essential features of the experiment when 
the external field H does not exceed the coercive field of 
the free layer. As a function of current density J, there 
is hysteretic switching between parallel and anti-parallel 
orientations of the free layer and the fixed layer. In the 
experiment, the critical currents for P^AP and AP— >P 
switching have opposite sign but are approximately equal 
in magnitude . This is not a feature of our T = 3 K phase 
diagram, but it is much more nearly true in our simu- 
lation at 3000 K, where thermal fiuctuations are large 
enough to mimic the effect of the (slow) current sweep 
rate used in the experiment. The scale we calculate for 
H is about half as large as seen in the experiment. 

High-Field Behavior: The macrospin model correctly 
models noise-driven, low amplitude, in-plane precession 
when H is larger than the coercive field and the cur- 
rent density is low. The existence of large-amplitude, in- 
plane and out-of-plane precession at higher J agrees qual- 
itatively with observation, but the precession frequency 
function oj{J) is not monotonically decreasing as found in 
experiment. The simulation predicts peaks in dV/dl as- 
sociated with telegraph-noise switching between two (or 
more) states of magnetization. These peaks are present 
in the experiment, as is two-state telegraph noise. We 
find an OPP/AP bistable phase field that is occupied in 
the experiment by a microwave quite "W" -phase. Micro- 
magnetic simulations by others suggest that that vortex 
creation and annihilation occurs in this phase field. 

Other Effects: The phase diagram of the sine-torque 
model differs mostly in detail with the phase diagram of 
the minimal model. An exception is the presence of a 
high-field, high current "fixed" phase for the sine torque 
where the macrospin freezes into a fixed angle with re- 
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spect to the fixed layer. This phase does not seem to 
occur in micromagnetic simulations. Another exception 
is the relative persistence of the OPP phase when we scan 
from high values of H to low values of H . This behav- 
ior is also absent from the micromagnetics simulations. 
Spin-pumping and current-induced magnetic fields do not 
change the phase diagram in any significant way. Large 
topologically changes do occur if we allow the Gilbert 
damping parameter to change with angle, but we have 
no guidance for the form a{9) should take. This seems 
like a fruitful direction for future research. 
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APPENDIX: ENERGY EXPRESSION 

The energy of the free layer includes a Zeeman energy 
Ez from the external field i?, a magnetostatic shape 
anisotropy energy Eg^ and a surface anisotropy energy 
(parameterized by K^) that vanishes in the limit that 
the free layer thickness c? — > oo. 
The Zeeman energy is 



Ez = VnoM.xa ■ H. 



The shape anisotropy energy is 



(16) 



Es - yi^oMsrii • Hd - yi/xoA/^m ■ AA • m, (17) 



where = M^M ■ rh and M are the demagnetization 
field and demagnetization tensor, respectively. Referring 
to Fig. the total energy is 



E = i/ioM.Viicos^e + Msin^eisin^t^ 
+N sin^ B cos^ 0] - ^loM^VH cos 6 
; — sm COS 



(18) 



where L,AI,N are the demagnetization factors for the 
£, y, X directions. These terms can be combined to give 



2E 



/iL cos^ 9 + Hm sin^ 9 sin^ 



-f/ijvsin^6'cos^(?!>-2ftcos6', (19) 
where h — H/Mg, and 



Hi. AK 
/..^L-^, /.M = M, h.=N-^^. (20) 

If we model the thin free layer as a very flat ellipsoid with 
semi-axis a > 6 ^ c, Eqs. (2.23-25) in Ref . 113 give 



Qj 6 



iV = 1 



a e2(l-e2)i/2 
cE 



a(l-e2)i/2^ 



(21) 
(22) 
(23) 



where K and E are com plete elliptic integrals with ar- 
gument e = a/I — If^ /a^. For the nominal geometry of 
Ref. Ull we have 2a = 130 nm, 26 = 70 nm, 2c = 3 nm, so 
L w 0.017, M « 0.035, N « 0.948. 
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